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(57) ABSTRACT 

A sensing device includes an elongated whisker element hav- 
ing a flexible cantilever region and a base region where a 
change in moment or curvature is generated by bending of the 
cantilever region when it contacts an object. One or more 
sensor elements cooperatively associated with the whisker 
element provide one or more output signals that is/are repre- 
sentative of two orthogonal components of change in moment 
or curvature at the whisker base region to permit determina- 
tion of object distance, fluid velocity profile, or object contour 
(shape) with accounting for lateral slip of the whisker element 
and frictional characteristics of the object. Multiple sensing 
devices can be arranged in arrays in a maimer to sense object 
contour without or with adjustment for lateral slip. 
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SENSING DEVICE WITH WHISKER 
ELEMENTS 

This application is a division of U.S. application Ser. No. 
11/906,751 filed Oct. 3, 2007, now U.S. Pat. No. 7,774,951, 5 
claims benefits and priority of provisional application Ser. 

No. 60/849,400 filed Oct. 4, 2006, the disclosure of which is 
incorporated herein by reference. 

CONTRACTUAL ORIGIN OF THE INVENTION 10 

This invention was made with government support under 
Grant No. IOB0446391 awarded by the National Science 
Foundation and Grant No. 1261628/NAS-03001 awarded by 
the National Aeronautics and Space Administration (NASA) 1 5 
through the Jet Propulsion Laboratory. The government has 
certain rights in the invention. 

FIELD OF THE INVENTION 

20 

The invention relates to sensing methods and devices that 
include one or more whisker sensing elements for detecting, 
sensing, measuring or quantifying object distances, contours 
and fluid flow profiles. 

25 

BACKGROUND OF THE INVENTION 

Facial and body whiskers serve as a remarkably sensitive 
source of tactile information for many species of mammals. 

For example, harbor seals can use their facial whiskers to 30 
track hydrodynamic trails, an ability thought to aid long- 
distance prey tracking (Dehnhardt et al. 2001). As another 
example, behavioral experiments on rats have shown that 
individual whiskers provide the fine-grained distance dis- 
crimination sensitivity needed to sense aperture width (Krupa 35 
et al. 2001). Increased study of the biological function of 
whiskers has been paralleled by an increasing interest in 
constructing their robotic counterparts (see references). 
Robotic whisker arrays could be used in a wide variety of 
applications, ranging from fine sensing tasks that push the 40 
limits of tactile discrimination and hence require engineers to 
construct systems that rival animals’ exquisite tactile sensi- 
tivity, to obstacle-avoidance tasks that require only far coarser 
sensing capabilities. 

The simplest types of artificial whiskers are those which 45 
are used as binary contact detection sensors. These have been 
used on several successful toys (e.g. “BioBugs” made by 
Wow Wee/Hasbro). Other researchers have explored the use 
of whiskers for wall-following and to characterize surface 
texture and surface defects. Most recently, a whisker sensor 50 
was designed for precise three-dimensional measurement of 
heart position in robot assisted beating heart surgery. 

An important feature of biological whiskers is their ability 
to extract three-dimensional (3D) features, either of solid 
objects or of fluid flows. However, few studies have investi- 55 
gated how this capability might be replicated in an artificial 
whisker array. One of the most successful approaches towards 
3D feature extraction was taken by Kaneko et al., in IEEE, 
Trans. Robotic Autom, 14, 278-29 (1998). These authors 
employ a method in which a flexible beam is rotated a small 60 
amount (“tapped”) against an object while measuring bend- 
ing moment at the base of the beam (whisker) to determine 
contact distance based on the rotational compliance. How- 
ever, this method has serious drawbacks: it requires multiple 
adjustments of actuator orientation to keep the beam oriented 65 
perpendicular to the object, to avoid lateral slip, and also 
requires multiple rotations for each radial distance extraction. 
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These are not only awkward in practice, but also infeasible 
when arrays of multiple whiskers are employed to contact the 
object. 

There are two interrelated problems that have as yet pre- 
vented artificial whiskers from being used in large, highly 
parallel, actuated arrays to sense object features. The first 
problem, as mentioned above, is lateral slip, in which the 
whisker slides out of its primary plane of rotation. The second 
problem is closely related to lateral slip, and involves the need 
to obtain a reasonable estimate of the coefficient of friction. 
The invention described herein is advantageous to overcome 
these problems and to provide and the ability to quantify and 
passively accommodate for lateral slip of a moving whisker in 
the presence of friction. 

SUMMARY OF THE INVENTION 

An embodiment of the invention provides a sensing 
method that involves relatively moving an elongated whisker 
element having a deflectable cantilever region, and an object 
so that the cantilever region bends as a result of contact with 
the object, and determining two orthogonal components of a 
change in a whisker bending parameter such as, for example, 
a change in moment exerted at a base region of the whisker 
element, or a change in curvature at the base region of the 
whisker element, as a result of the bending wherein the 
method further includes estimation of object distance (con- 
tact point distance) by accounting for lateral slip of the whis- 
ker element. The estimate of contact point distance is adjusted 
by a factor related to the object surface slope relative to the 
plane of movement of the whisker, and frictional properties of 
the object/whisker contact surface. An array of multiple whis- 
ker elements and the object can be relatively moved in prac- 
tice of the invention. 

A sensing apparatus pursuant to an embodiment of the 
invention includes a sensing device having an elongated 
whisker element with a flexible cantilever region and a base 
region where a moment is generated by deflection of the 
cantilever region, and one or more sensor elements coopera- 
tively associated with the whisker element to provide one or 
more output signals that is/are representative of two orthogo- 
nal components of a change in a whisker bending parameter 
(e.g. change in moment or change in curvature) exerted at the 
whisker base region. The output(s) of the one or more sensing 
elements is/are input to a computer device having an algo- 
rithm which uses the values of the two measured orthogonal 
components of change in moment or change in curvature to 
estimate object distance by accounting for lateral slip of the 
whisker element along the object without any need to physi- 
cally reorient the whisker element or the object or to sense 
force. Multiple sensing devices can be arranged in arrays in a 
manner to sense distance or velocity profiles or object shapes. 

In a particular embodiment, the base region of the whisker 
element is mounted on a mounting body. One or more sensor 
elements is/are mounted on the mounting body in manner to 
sense change in moment or change in curvature at the base 
region from forces acting on the flexible region of the whisker 
and transmitted to the mounting body. Preferably, multiple 
sensor elements are disposed on the mounting base to provide 
one or more output signals that are representative of two 
orthogonal components of change in whisker moment or 
whisker curvature exerted at the whisker base region from 
bending of the whisker. 

Certain apparatus embodiments and method embodiments 
of the invention are practiced in a so-called tapping mode 
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where the whisker element is moved (rotated and/or trans- 
lated) through a small displacement against an object to deter- 
mine initial contact distance. 

Certain further system embodiments and method embodi- 
ments of the invention are practiced to detect fluid flow pro- 
files. 

Other features and advantages of the invention will become 
apparent form the following description taken with the fol- 
lowing drawings. 

DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic view of sensing device having a 
whisker element rotated by motor against an object (peg). 

FIG. 2 a is a schematic view of the whisker element rotated 
against an object and illustrating longitudinal slip. 

FIG. 2b is a schematic three dimensional view of the sens- 
ing device illustrating the plane of rotation (the x-y plane), 
lateral slip, and the sensing plane that geometrically describes 
lateral slip. 

FIG. 2c is a schematic of a sensing plane used in the 
calculation set forth in the APPENDIX. 

FIG. 3 is a schematic perspective view of a sensing device 
pursuant to an illustrative embodiment of the invention. 

FIG. 4 is a perspective view of the sensing apparatus used 
for model validation experiment of Example 1 . FIG. 4a is a 
simplified view of the apparatus showing the angles involved. 

FIG. 5 is a perspective view of the sensing apparatus used 
for model validation experiment of Example 2. 

FIG. 6 is a complete splined surface and associated contact 
points for Example 2. 

FIG. 7 is a front view of the sensing apparatus used for 
sensing fluid flow profile of Example 3. 

FIG. 8 is a graph of fluid flow profile of air velocity versus 
whisker (strip) height in the array of Example 3 where mea- 
surements were made with a compressed air source at two 
different distances in front of the array and where triangle 
data points are obtained from whiskers and circle data points 
are obtained from Pitot tube. 

FIG. 9 is a schematic side view of three bundles of whis- 
kers, each containing three whiskers of different length. 

DESCRIPTION OF THE INVENTION 

Classical elasticity theory was used to model the whisker 
element W as a cylindrical beam bending against a rigid, 
motionless object (shown as peg P) and to derive the mono- 
tonic relationship between radial contact distance d, the push- 
ing angle 0 and moment M (torque) at the whisker base (see 
FIG. 1). Through either rotation or translation (or a combina- 
tion of these motions) of the whisker at or near the base region 
thereof (where the moment sensor resides) against an object, 
the initial point that the whisker touches along the object 
during the whisker’s movement can be computed as 
described in the APPENDIX. The derivation can alternately 
be conducted in terms of force instead of moment, and/or 
translation instead of rotation. Although the invention is 
described in detail immediately below with respect determin- 
ing a change in moment at the base region of the whisker 
element as a result of bending, this is offered for purposes of 
illustration and not limitation since the invention is not so 
limited and can be practiced by determining a change in 
curvature or other deflection of the whisker element at the 
base region (or other region) as explained later below. The 
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invention envisions moving the whisker element, or the 
object, or both to provide relative motion there between. 

From the APPENDIX, equation 2 shows that a motor rotat- 
ing a cylindrical beam against a point object will experience 
5 a torque proportional to the amount it is rotated (for small 
values of 0) ifthereisno lateral slip. In other words, the beam 
acts as a torsional spring, with an effective spring constant 
inversely proportional to object distance d. Kaneko et al. in 
IEEE, Trans. Robotic Autom, 14, 278-29 (1998) reported 
10 such a proportionality. As described above, however, the 
reported method of Kaneko et al. is disadvantageous in that 
this method requires adjustment of actuator orientation over 
multiple whisks to orient the whisker oriented perpendicular 
15 to the object being sensed and avoid lateral slip. 

An illustrative embodiment of the present invention pro- 
vides a sensing device that can provide appropriate output to 
an algorithm of a computer device that accounts for lateral 
slip of the whisker W without the need to adjust actuator 
20 orientation to keep the whisker W within the plane of its 
rotation. 

Referring to FIG. 2a, longitudinal slip is illustrated for a 
whisker element W rotated by an electric or other motor 
shown where angular motion of the motor shaft is detected by 
25 angle sensor, and moment at the base region of the whisker W 
is sensed by moment (torque) sensor. Longitudinal slip is a 
relatively minor complication that arises when the point of 
object contact is not a sharp edge, but instead has a finite 
curvature in the plane of rotation (longitudinal curvature) as 
11 1 depicted in FIG. 2a. The result is a decrease in the resulting 
rotational compliance and underestimation of contact dis- 
tance. However, the decreased compliance will be negligibly 
small unless the object surface has extremely low curvature or 
35 contact occurred very close to the whisker base. A much more 
significant complication arises when the object surface is 
slanted (inclined) relative to the plane of whisker rotation, as 
depicted in FIG. 2b. That is, the object surface is not perpen- 
dicular to the plane of rotation. In such cases, a phenomenon 
40 called lateral slip occurs where the point of contact slides 
significantly along the periphery of the object, the whisker 
bends outside the plane of rotation and the resulting compli- 
ance is much greater than would have occurred without slip, 
causing equation 2 of the APPENDIX to overestimate d. 

45 As described above, the reported method of Kaneko et al. 
in IEEE, Trans. Robotic Autom, 1 4, 278-29 (1998) adjusts the 
actuator orientation, based on measurement of out-of-plane 
moment M v , to keep the plane of rotation oriented perpen- 
dicular to the object being sensed, avoiding lateral slip. This 
50 permits exploitation of the direct proportionality rule that 
exists when the whisker deflects within the plane of rotation. 
However, this method is disadvantageous due to the need to 
actively reorient the whisker multiple times. 

Accounting for the lateral slip of the whisker W requires 
the ability of the torque sensor to measure the bending 
moment M r out of the plane of rotation. However, instead of 
actively adjusting the plane of rotation to avoid lateral slip, an 
estimate of the friction coefficient \i est that acts between the 
60 whisker element and the object is used to estimate lateral slip 
and hence the object distance d. 

Pursuant to an embodiment of the invention, the lateral slip 
is sensed and its presence is passively accounted for by mea- 
suring out-of-plane moment M F , and utilizing a reasonable 
65 estimate of the friction cone angle between the whisker ele- 
ment and the object surface. In practice of this illustrative 
embodiment of the invention, equation 2 of the APPENDIX is 
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modified to incorporate moment M v and friction cone angle 
estimate a est , yielding equation 15: 


d -/ : ' ' J 

A M z + Mytan/3 est 

where k = 3/;7 and = a est + taiT l (M y / M z ) 

(see APPENDIX). 

10 

E is elastic modulus of the whisker element, and I is the 
moment of inertia of the whisker element. The estimate of 
object distance thereby accounts for lateral slip in equation 1 5 
in this embodiment by measuring and accounting for deflec- j - 
tion of the whisker element outside of the plane of movement. 
The results [estimated object distance(s) and/or contours] 
determined in practice of the invention are stored, displayed, 
or conveyed to or available to a user (e.g. see FIGS. 6 and 8). 

The need to adjust the plane of rotation is thus effectively 20 
eliminated by practice of the invention as long as a reasonable 
estimate can be made for the friction coefficient between the 
whisker and the surface. As explained in the APPENDIX, 
inaccurate estimate of friction cone angle between the whis- 
ker element and the object surface can lead to distance extrac- 25 
tion inaccuracy. However, so long as a es , (friction cone angle) 
falls within a moderate range of true friction cone angle, the 
distance extraction will be accurate even for steep surface 
slopes. For example, if a ejt =10°, distance extraction will be 
accurate to within 5% so long as the true friction cone angle 30 
a ranges between 8.2° and 1 1.7° and surface slope (I is less 
than or equal to 60°, or to within 10% accuracy so long as 
friction cone angle a ranges between 7.8° and 12.3° and 
surface slope (I is less than or equal to 70°. 

With respect to lateral curvature of the object surface, the 35 
model of the APPENDIX assumes that the surface being 
contacted is flat; i.e. has a lateral curvature of zero. The effect 
of curvature on distance extraction must be considered in the 
context of three additional variables: 0, a and fl 0 (surface 
slope at the initial contact point). When so considered, the 40 
consequence of lateral curvature of the surface being con- 
tacted is relatively mild. For example, given a significantly 
curved surface with a normalized radius of curvature of the 
surface, r„=l (i.e., lateral radius of curvature equals radial 
contact distance d), and a reasonable pushing angle of 0=3°, 45 
a surface slope |3 0 of 63° will result in a 5% distance extraction 
error and a surface slope (3 0 of 74° will result in a 10% 
distance extraction error. Increasing 0 to 5° would decrease 
the maximum surface slope |3 0 to 52° and 65°, respectively. 

Referring to FIG. 3, a sensing device is shown pursuant to 50 
an illustrative embodiment of the invention that can supply 
appropriate output to an algorithm of a computer device that 
incorporates equation 15 to accomplish radial distance 
extraction. The sensing device includes an elongated whisker 
element W having a flexible cantilever region CR with a free 55 
end FE and a base region BR where a moment is generated by 
deflection (bending) of the cantilever region as it contacts the 
object. One or more sensor elements G are cooperatively 
associated with the whisker element W to provide one or more 
output signals that is/are representative of two orthogonal 60 
components of moment, M v and M_ exerted at the whisker 
base region. 

In particular, in FIG. 3, an artificial cylindrical whisker 
element W has a base region BR that is fixedly mounted on a 
mounting body B to provide a fixed center of rotation at a 65 
fixed speed of rotation. The one or more sensor elements G 
are mounted on the orthogonal faces of the mounting body B 
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in maimer to sense two-axis moments M_ and M v at the base 
region BR from forces acting on the flexible cantilever region 
CR of the whisker and transmitted to the mounting body. 

In an illustrative embodiment of the invention offered for 
purposes of illustration and not limitation, the sensing device 
can comprise a straight superelastic Nitinol metallic wire (1 
nun diameter, 1 0 cm length), with a small set screw SI (FIGS. 
4 and 5 ) attached at the base region BR of the Nitinol wire by 
adhesive (e.g. cyanoacrylate-superglue) or any other attach- 
ment mechanism. The set screw SI allows different size and 
shape Nitinol wires to be easily interchanged within the 
mounting body B. The mounting body B comprised a small 
aluminum block (4x4x8 nun) with the center tapped to 
receive the outer threaded region of the set screw SI. Each of 
the four flat faces of the mounting body was fitted with a 
conventional strain gauge G, allowing independent measure- 
ment of M„ and M v , and the mounting body B itself is attached 
to a larger set screw S2 to allow easy connection to an array 
shown in FIG. 5 for purposes of illustration and not limitation. 
In FIGS. 4 and 5, the array design comprised a vertically 
oriented aluminum plate P with a series of threaded holes H 
that allows any given number of sensing devices having set 
screws S2 to be mounted and configured in a single colunm. 
,\n aspect of this particular embodiment of the invention, 
although it is not necessary for all embodiments, is that it 
positions the base of the mounting body (and the fixed tips of 
their cantilever regions) at the center of rotation. This simpli- 
fies distance extraction because the analysis exactly conforms 
to the model described in the APPENDIX; however, the equa- 
tions can easily and obviously be modified to accommodate 
for rotations and/or translations when the whisker base is not 
at the cent er of rotation. 

The invention envisions using strain gages or other alter- 
native two-axis deflection and/or two-axis moment and/or 
two-axis curvature sensor elements that include, but are not 
limited to, light emitting diodes/receivers, fiber optics, optical 
emitter/detector pairs (described in reference Hirose et al. 
1990), magnetic sensors including Hall-effect sensors (de- 
scribed in reference Hipp et al. 2006), capacitive coupling 
(described in reference Fend et al. 2003), linear differential 
transformer (LVTD) sensors, video cameras (described in 
Kaneko et al. 2001), electric dial gages (described in refer- 
ence Kaneko and Tsuji 2000), and piezoelectric sensors. The 
sensor elements preferably are selected to determine two 
orthogonal components of a change in moment, curvature or 
other deflection of the whisker element as a result of its 
bending. For example, the moment of the whisker element at 
the base region is directly proportional to the change in cur- 
vature of the whisker element at the base region (the whisker 
element, including its base region, need not be straight while 
at rest). Specifically, Ak=M(s)/EI where Ak is change in 
curvature, s is the arc length, M is moment, E is elastic 
modulus, and I is moment of inertia of the whisker element 
wherein E and I may also vary with s. The above mentioned 
sensor elements can be used in practice of the invention to 
determine two orthogonal components of a change in curva- 
ture (or moment) of the whisker element at the base region as 
a result of bending of the whisker when it contacts an object. 

The output of the sensor elements is representative of the 
two orthogonal components of change in curvature (or 
moment) of the whi sker element for use as described above as 
input to the computer to determine object distance with 
accounting for lateral slip. 

In practicing the method described above, consistently reli- 
able object distance extraction in the presence of significant 
surface slopes includes a reasonable estimate of the friction 
coefficient between the whisker element and the obj ect. There 
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are three possible methods by which an estimate of friction 
mi gilt be obtained, as follows: 

1 . If the orientation of the entire whisker array can be tilted 
along the x-axis of any particular whisker, whisks would 
be performed against the object at a variety of orienta- 
tions, effectively adjusting the surface slope. The range 
of slopes for which (j>=0° would be equal to 2 a (the 
accuracy depending on how gradually the orientation 
was adjusted). This essentially amounts to an in-field 
test of the object to determine the friction coefficient. 

2. Have the array explore the object, and repeat the feature 
extraction algorithm (the transformation of contact 
points into 3-D object shape, e.g. splining) using a range 
of friction coefficients. The one that leads to the most 
consistent contact point positions and thus the smoothest 
surface is kept. 

3. Have the array whisk against the object at a variety of 
array orientations. All the contact points sampled are 
given a confidence measure related to <(>, such that points 
with large 4> are ignored or have less influence on the 
feature extraction algorithm. 

Note that methods 1 and 2 assume the friction coefficient is 
consistent over the entire object, while method 3 does not. 

One final issue involves contact occurring at the whisker 
free tip or end. Such cases may result in increased rotational 
compliance, causing equation 1 5 to estimate d r to be greater 
than the whisker length. The solution is to simply give d T a 
threshold equal to the whisker length. It is also worth noting 
that regions of an object that are concave with respect to the 
plane of whisker rotation can only be sampled through tip 
contact with a straight whisker. Tip contact is thus actually 
desirable for two reasons: it potentially helps provide a very 
precise measurement of contact distance (since the whisker 
length is fixed), and allows concave regions of an object to be 
sampled. 

Example 1 

With respect to FIGS. 4 and 4a, an experiment involved 
rotating the whisker W against a slender stainless steel bar at 
a variety of surface slopes (see FIG. 2b): from surface 
slope=0° (vertical) to 75° (nearly horizontal), in 5° incre- 
ments, and at two different speeds: =10 deg/s and 90 deg/s. 
Surface slope was determined using the protractor shown. 
Signals from the strain gages at the base were first calibrated 
from voltage to moment. The M v component was calibrated 
by rotating the whisker against the bar with (3=0° at a radial 
distance of d x =5 cm. Ten whisks were performed at each of 
the two velocities. The location at which the whisker gently 
touched the bar was defined as 0=0°. Whisker trajectories 
were chosen to ensure a smooth acceleration to the maximum 
velocity. For 0=10 deg/s the whisker was rotated through a 
range of -3° to 12°, while for 0=90 deg/s, an expanded range 
of-15°to20° was used. Hie mounting body was then rotated 
90° and the calibration process repeated for M_. 

All data were filtered at 800 Hz, sampled at 2000 Hz and 
passed through a zero-phase digital filter with a cutoff fre- 
quency of either 5 Hz (for 0=10 deg/s) or 45 Hz (for 0=90 
deg/s). The strain gages each were incorporated into a Wheat- 
stone half-bridge. All data analysis was performed using 
MATLAB v7.0 software on an Intel Pentium IV computer 
having equation 15 programmed as part of an operating algo- 
rithm. Equation 2 of the APPENDIX with a push angle 0=3° 
was used to arrive calibration factors between voltage and 
moments M v and M_. 

Experimental data were gathered in the same way as the 
calibration process, but for varying surface slopes. For each 
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surface slope, the plane of the bar was perpendicular to the 
orientation of the whisker at 0=0°. An excellent fit between 
model (equations 9, 10 and 11) and experiment was obtained 
for whi sking against the steel rod at 1 0 deg/ s, with the result- 
5 ing estimation of friction cone angle being a ei< =10.5°. At 90 
deg/s, the data again conforms well to the model, the fit 
yielding a CJ( =12.3°, but the standard deviations are signifi- 
cantly higher. A rough surface obtained by providing 1 80 grit 
sandpaper stuck on the stainless bar introduced a problem 
to regarding friction known as stick-slip, in which the whisker 
would irregularly switch between periods of no slip, smooth 
slip and sudden, abrupt slip. This led to moderately inconsis- 
tent results for 0=10 deg/s. However, stick-slip seemed to 
have a less significant impact for 0=90 deg/s. Hie fits to 0=10 
15 and 90 deg/s. yielded a erf =23.3° and 18.7°, respectively. 

The above-described embodiment of the invention pro- 
vides a so-called tapping method for implementing radial 
distance detection with an artificial whisker W in the presence 
of significant lateral slip (lateral deflection in the z-direction) 
20 and surface-friction. By modeling the whisker as a rotating 
(or translating) cantilever beam, it is possible to interpret the 
sensed moment to infer information about the object being 
touched. Using tapping techniques, the whisker need only 
rotate (whisk) a small (approximately 2 to 5°) angle past 
25 initial object contact, so that the radial distance of contact can 
be inferred. The method requires only that the whisker be 
equipped with a two-axis sensor at the base that can measure 
moment and/or curvature, and works best when a reasonable 
estimation for the object’s friction coefficient is available. An 
30 important advantage of this method is that it does not require 
that each whisker can independently adjust its plane of rota- 
tion, and is thus very amenable to implementation on a large- 
scale array. 

35 Example 2 

This example involves a whisker array for quickly and 
consistently extracting complex object features using only a 
small number of whisks. To investigate the shape-extraction 
40 capabilities, a 4x1 vertical array of robotic whiskers W was 
constructed as shown in FIG. 5. [space inserted here] Each 
Nitinol wire cantilever region CR emerged from a set-screw 
SI (visible in the figure), which was inserted into a rectangu- 
lar mounting body (B), also called a “follicle.” The follicle 
45 was fitted with four strain gages on the faces (not shown so see 
FIG. 3). In this embodiment, the follicle constitutes the whis- 
ker base region used to measure the two orthogonal compo- 
nents of moment as described above. The base of each follicle 
was threaded to match the threads in the holes on the 4x1 
50 array, so that each follicle could easily be inserted into the 
array. 

The vertical array was tested by tapping (“whisking”) it 
against a small sculpted head H, shown in FIG. 5 — a complex 
object selected specifically for its intricate concavities and 
55 convexities. In particular, the plate P was rotated by an AC 
motor about a vertical rotational axis that included the fixed 
end regions of the whisker elements W. Whisks were per- 
formed with the array positioned at evenly spaced hei gilts and 
angles relative to the sculpted head so as to cover the head’s 
60 entire left side. Each whisk provided four radial contact 
points in 3D space (one for each whisker). Splining the con- 
tact points and mirror-imaging the left side to the right then 
resulted in a faithful extraction of the original shape of the 
sculpture, FIG. 6. 

65 In particular, calibration between voltage and moment M v 
for the artificial whisker array was performed by sweeping 
each whisker five times against a peg placed at a single dis- 
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Example 3 
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tance of 2 cm. In this particular example, a low whisking 
velocity (10 deg/sec) minimized inertial effects and maxi- 
mized the resolution of contact angle detection. These param- 
eters could be changed to be more suitable for any particular 
application. This same procedure was then done with the 5 
whisker mounting body rotated 90° to calibrate between volt- 
age and moment M z . 

To gather experimental data, the sculpted head was fixed at 
the center of a cylindrical coordinate system (r, 0, z) while the 
array whisked at several different positions. Note that 0 of this 10 
general coordinate system is not the same as the angle 0 
through which the whisker deflects. The array was positioned 
at regular intervals of height z and angle 0. The choice of 
radial distance from the motor shaft to the head was neces- 15 
sarily dictated by the requirement that the whiskers make 
contact with all regions of the left side of the head (Table 1 ). 

TABLE 1 

The array was positioned in regular intervals of height z and angle 0, ^ 

while distance r was manually chosen to ensure whisker contact 
with all regions of the face. 


z (cm) 0 (degrees) r (cm) 


0.000-2.000, 0.125 cm increments 0-60°, 10° increments 5.00 25 

2.125-3.250, 0.125 cm increments " 4.75 

3.375- 4.250, 0.125 cm increments " 5.50 

4.375- 6.000, 0.125 cm increments " 5.25 


A single whisk was perforated at each position. Analog 30 
signals from each whisker base were filtered at 160 Hz, 
sampled at 500 Hz and then passed through a zero-phase 
digital filter (period =1 /2 second). All data analysis was per- 
formed using MATLAB v7.0 software on Intel Pentium IV 
computer (CPU) having equation 1 7 programmed as an oper- 35 
ating algorithm (i.e., a esr was assumed to be zero). 

The location of each contact point was then computed 
relative to the array. The angular component 0 was taken to be 
the angle at which the whisker first made contact with the 
object (when dIMI/dt crossed a small threshold), and the 40 
radial d component was found using equation 17 of the 
APPENDIX. 

A total of 343 whisks were perforated. In cases where a 
whisker did not make contact, or the data analysis showed that 
the whisker hit on the right side of the face, the data were 45 
automatically removed. All other contact points were mirror- 
imaged to the right side of the head. The longer whiskers 
captured the broad convex features, while the shorter whis- 
kers explored the finer concave regions. Equation 1 7 was also 
used to determine when contact occurred at a whisker tip . Tip 50 
contact always resulted in an estimate of d very close to or 
greater than the whisker length. In the case that estimated d 
was greater than the whisker length, d was automatically set 
equal to the length. Thus, equation 1 7 provided the informa- 
tion that contact did indeed occur at the whisker tip and not 55 
mid-length. This would not have been possible had the sens- 
ing mechanism been, for example, a simple binary contact- 
switch. 

Note that regions of the sculpted head that are concave in 
the plane of rotation (for any given height) are only reachable 60 
through tip contact, which accounts for 25% of the splined 
surface. A total of 1036 contact points on the left side were 
collected, with 438 points determined through equation 1 7 to 
occur at or near the tip (42% of all points). The points were 
converted to head-centered C artesian coordinates to simplify 65 
splitting, mirror-imaged to the right side, and plotted along 
with the spline to create FIG. 6. 


To examine the ability of whiskers to characterize fluid 
flows, two opposing 4x1 whisker arrays shown in FIG. 7 were 
used. To estimate the shape of a variable fluid flow profile in 
two dimensions, two artificial whisker arrays were arranged 
in a plane. The sensing element was reduced to one dimension 
(i.e., M_ only, although M v can also be used in general) and the 
whiskers took the form of highly-flexible plastic strips to 
increase the surface area exposed to the flow and maximize 
bending. The plastic strips were Flex Sensors, a commercial- 
off-the-shelf sensing technology that can measure bending 
along the plastic strip. However, in this embodiment, all bend- 
sensitive elements along the length of the Flex Sensor were 
shorted out except the one sensor element G closest to the 
base. The flexible strips are glued to the vertical stainless steel 
support rods R having a length of 1 9 cm. The dimensions of 
the each flexible strip are 0.5 cm width by 1 1 cm length, the 
strips being available from Jameco Electronics. Compressed 
air is directed from a compressed air source towards the 
center of the plane formed by the flexible strips. 

The sensor element of each flexible strip was calibrated 
over a range of deflections to account for any nonlinear sens- 
ing properties, giving a calibration factor from voltage to 
moment. A compressed air source was centered at two dis- 
tances ( 1 5 cm and 40 cm) in front of the array, and data were 
averaged over five seconds for each whisker. Averaged analog 
voltage signals from the strain gages at each whisker base 
were directly converted to moment using the calibration fac- 
tor. All data analysis was perfonned using MATLAB v7.0 
software on an Intel Pentium IV computer programmed with 
the operating algorithm described in the next paragraph 

First, the true value of the air speed was determined by 
placing a Pitot tube with an attached manometer at the height 
of each whisker. Next, the square root of moment (moment 
being proportional to distributed load) was calibrated to air 
speed to characterize the flow, consistent with the Bernoulli 
equation for incompressible flow. Ten trials were taken for 
both the experimental (whisker) and control (Pitot tube) data. 

When the stream of air was directed toward the center of 
the array, each whisker deflected by a magnitude dependent 
on its distributed load, permitting accurate characterization of 
the shape of the fluid stream. As shown in FIG. 8, velocity 
estimates using the (square root of) moment from the base of 
each whisker matched up will with the Pitot tube data. More 
complex whisker configurations would allow flow patterns to 
be described in three dimensions, potentially enabling wake 
following (remote tracking) of moving underwater objects. 

It is noted that this technique does not give information 
about the load as a function of position along the whisker 
length. Instead, each whisker in the array gives a single esti- 
mate of flow velocity at a given height by integrating the 
distributed load over its length. Because the eight whiskers 
are at different heights, the flow profile can be characterized 
in the z-direction. Using whiskers of different lengths and in 
different array configurations could allow three-dimensional 
characterization of the flow pattern. 

Underwater applications will benefit from the ability to 
sense fluid flow profiles. Again using a traditional cantilever 
beam model, but also arranging bundles of several whiskers 
of different length in parallel and with their base regions 
closely aligned as in FIG. 9, it is possible to estimate the ID 
(one dimensional) velocity profile of the fluid flow along the 
bundle’s axis. Arranging several such bundles of whiskers in 
a plane, it is possible to estimate the velocity profile of the 
fluid within the cross-section of the plane to arbitrary accu- 
racy (depending on the number of whiskers and their spac- 
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ing). Furthermore, arranging several such cross-sections of 
whiskers in a row can allow for description of the velocity 
profile within a volume of fluid. Practical considerations of 
this method may involve whiskers colliding with one another, 
and specifically for the three-dimensional case, whiskers 
planes near the front of the flow influencing the flow profile 
before it reaches whisker planes that are towards the back of 
the flow. These results have potential applications to both 
land-based robots and underwater autonomous vehicles, in 
which robust tactile perception could enable a wide range of 
behavioral capabilities. 

APPENDIX 

Calculations 

The whisker is modeled as a straight, cylindrical, flexible 
beam rotating with a fixed center of rotation and at a constant 
velocity. At some point, the whisker comes into contact with 
an object, at which time the task of finding the distance from 
the whisker base to contact point arises. We assume the object 
is rigid and that the point of contact is fixed, discrete, and 
exists somewhere along the length of the whisker (as opposed 
to the tip, a condition which is discussed in the Description Of 
The Invention Section). In the case that lateral slip is negli- 
gible, the configuration of the whisker can be described in 
2-dimensions as shown in FIG. 1. 

Because only a small deflection is needed to sense the 
rotational compliance, classical linear beam theory can be 
used. The whisker rotates by a small pushing angle 0 against 
an object at radial distance d^, imposing a deflection d v . The 
resulting bending moment m z at the base of the whisker is 
measured by a torque sensor. Under these conditions, Kaneko 
et al. showed that the rotational compliance C e is directly 
proportional to object distance, and can be expressed as: 

<k=kCe (1) 

where k=3EI, E is the elastic modulus of the whisker and I is 
the area moment of inertia of the whisker (I=jc/ 4- 4 , where r is 
the radius of the beam’s cross-section). Rotational compli- 
ance is defined as the ratio of angular displacement to bending 
moment at the base, allowing us to write 



When the orientation of the object is not perpendicular to 
the plane of rotation, lateral slip may occur, in which case the 
contact point will drag along the periphery of the object, as 
shown in FIG. 2b. The resulting rotational compliance will be 
greater than when lateral slip does not occur, causing equation 
2 to overestimate d T . There are three specific conditions at the 
region of contact that affect if and how lateral slip occurs: 

1. the slope of the surface at the contact point within the 
sensing plane, 

2. the curvature of the surface within the sensing plane, and 

3. the frictional properties between the whisker and the 
object. 

In the results that follow, the slope of the surface will be 
directly accounted for. We know that lateral curvature has 
only a mildly detrimental impact on distance extraction in 
most cases as long as 0 is sufficiently small. Friction can 
present some difficulties, although there are ways to address 
them which are discussed. For now, we assume that the sur- 
face curvature is zero (i.e. is flat) and that traditional Coulomb 
friction exists between the whisker and the surface. Because 
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this is a quasi-static analysis, any possible discrepancies 
between static and kinetic friction coefficients are assumed to 
be negligible. 

To sense lateral slip, the torque sensor (strain gage(s)) at the 
5 base of the whisker must be equipped to measure out-of-plane 
bending of the whisker, i.e. bending moment M v . Understand- 
ing how to make use of M v to account for the lateral slip is 
aided by analysis of the sensing plane, illustrated in FIG. 2b. 

The sensing plane is an imaginary 2-D projection of the 
1 0 region where the whisker is touching the obj ect, parallel to the 
y-z plane and intersecting the contact point. In FIG. 2b, the 
dashed line indicates the initial contact point, the dotted line 
indicates the current contact point, and the dash-dotted line 
j - indicates where the whisker would intersect the sensing plane 
were it not obstructed by the surface — the so-called ‘virtual 
point’. A basic geometrical analysis of the relationships 
between variables in the sensing plane is sufficient to derive a 
modified form of equation 2 that accounts for the lateral slip. 
20 FIG. 2c defines the relevant variables projected on the 
sensing plane. It geometrically predicts the lateral slip a of the 
whisker along the object, depending on the local surface slope 
(I and the friction cone angle a. Slip will not occur when ag |3 
because in this case no movement of the contact point is 
25 required to keep the contact force f within or on the boundary 
of the friction cone, and the resulting slip angle <j> will be zero. 
Conversely, when a<|B, the contact point will slip just enough 
to enable static equilibrium, resulting in a contact force angle 
of (I-a. 

11 1 Summarizing, we have 

a>p.<p = 0 (3) 

a < ft, <p = ft - a 

which can be expressed more compactly as 

<|)=max(0,p-a). (4) 

40 

Notice that the overall linear deflection is now split into two 
components: d z and d . The distance between the virtual point 
and the initial contact point is simply the contact distance d,. 
multiplied by 0. We also can now give a concrete geometric 
45 definition to the lateral slip a as the distance the whisker slips 
along the object surface in the lateral plane: 

d x 9tan(p - (pcosft (5) 

5 Q tan0sin/? 

Using this model of lateral slip, we seek an equation analo- 
gous to equation 2 that will allow determination of the dis- 
tance d T based on known and measurable variables : E, 1, 0 , M, 
55 ... 4 

and M r . In the absence of lateral slip, it is clear that 0=tan(d y / 

d x ). For small angles, 0=d ) /d. c , which we plug into equation 2 
to obtain: 



Note that this equation holds when lateral slip occurs (FIG. 
65 2b) because it simply relates moment in the z-direction to 
deflection in the y-direction, neither of which are affected by 
independent deflection of the whisker in the z-direction (d z ). 
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By symmetry, the same linear scaling between m z and d 3 . also which can alternatively be written as 

holds for m v and d z : 


M 


y 



(7) 5 


d x = k cos d>, 

\M\ V 


( 17 ) 


Furthermore, from FIG. 2c, it is clear that 


where IMI=Vm z 2 +M v 2 . 
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Combining equations 6, 7 and 8 and solving for d r yields: 
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Unfortunately, there’s no way to reliably estimate the sur- 
face slope |3 in the presence of unknown friction. The most 4(j 
straightforward concession that can be made to arrive at a 
solution is to assume or estimate some nominal friction coef- 
ficient [i esr . First, we note that 


a est =tan 1 


(13) 


If we now assume <x es fik$ est , $\ est can easily be computed 
using equation 3: 

|3„,=<|)+a e „. (14) 

Replacing the unknown surface slope (5 from equation 1 2 

with P e ^y ields 
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Note that equation 1 5 is reliable even when a. es / > $ est (con- 
trary to equation 14’s assumption that a est =$ est ) because in 
this caseM v =0 and the tenn M v tan $ est becomes zero as well. 

One obvious situation to consider is that of frictionless 60 
contact, resulting in |U r =<t> and thus 
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described in detail above, those skilled in the art will appre- 
ciate that changes and modifications can be made therein 
within the scope of the invention as set forth in the appended 
claims. 15 

The invention claimed is: 

1. A sensing apparatus comprising: 

(a) a plurality of sensing devices, each said sensing device 
comprising: 

(i) an elongated whisker element having a deflectable 2 o 
cantilever region and a base region where a change in 
moment or curvature is generated by deflection of the 
cantilever region, and 

(ii) one or more sensor elements cooperatively associ- 
ated with the whisker element to provide one or more 2 5 
outputs that is/are representative of change in moment 

or curvature exerted at the base region; 

(b) a fluid flow source to direct fluid toward the plurality of 
sensing devices; and 

(c) a computer device to determine a fluid velocity profile 30 
using the one or more outputs of one or more sensing 
devices. 
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2. The sensing apparatus of claim 1, wherein the plurality 
of sensing devices are arranged in an array to determine fluid 
velocity profile along a dimension or axis. 

3. The sensing apparatus of claim 2, wherein each sensing 
device comprises whisker elements of different lengths 
arranged in a bundle. 

4. A sensing apparatus comprising: 

(a) a plurality of sensing devices, wherein each sensing 
device comprises: 

(i) an elongated whisker element having a deflectable 
cantilever region and a base region where a change in 
moment or curvature is generated by deflection of the 
cantilever region, and 

(ii) one or more sensor elements cooperatively associ- 
ated with the whisker element to provide one or more 
outputs that is/are representative of a change of two 
orthogonal components of moment, M and M_, 
exerted at the base region; 

(b) a fluid flow source to direct fluid toward the plurality of 
sensing devices; and 

(c) a computer device to determine a fluid velocity profile 
using the one or more outputs of the plurality of sensing 
devices. 

5. The sensing apparatus of claim 4, wherein the plurality 
of sensing devices are arranged in an array to determine fluid 
velocity profile along at least one dimension or axis. 

6 . The sensing apparatus of claim 5, wherein each sensing 
device comprises whisker elements of different lengths 
arranged in a bundle. 



